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Abstract 

The theory of Raman scattering by the electron-phonon coupled system 
in metals and heavily doped semiconductors is developed taking into account 
the Coulomb screening and the electron-phonon deformation interaction. The 
Boltzmann equation for carriers is applied. Phonon frequencies and optic 
coupling constants are renormalized due to interactions with carriers. The 
k— dependent semiclassical dielectric function is involved instead of the Lindhard- 
Mermin expression. The results of calculations are presented for various values 
of carrier concentration and electron-phonon coupling constant. 

PACS number(s): 63.20.Dj, 63.20.Kr, 71.38-k, 72.30.+q, 78.30.-j 



1 Introduction 

Recently there has been considerable interest in the effect of electron-phonon interac- 
tions on the optical-phonon dispersion. This interest is stimulated by contradictions 
between different approaches to the electron-phonon interaction. The strong phonon 
renormalizations were obtained first by Migdal pQ (see also |2]) within a consistent 
many-body approach based on the Frolich Hamiltonian. The extremely large dis- 
persion of optical phonons was predicted in work [3| using also the Frolich model. 
These results contradict to the Born-Oppenheimer (adiabatic) concept [3] according 
to which the phonon renormalizations should be small in terms of the nonadiabatic 



parameter ym/M, where m and M are the electron and ion masses, respectively 
(see, also [5]). Theoretical investigations ;6J of the sound velocity and acoustic at- 
tenuation in metals confirm the adiabatic concept. In a recent paper, Reizer [7j 
emphasized the importance of screening which should be taken into account. To 
our knowledge, the Coulomb screening effect on the LO phonons was firstly stud- 
ied in the work [HJ. In our work using the Boltzmann equation we have found 
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that the electron-phonon interaction results more essentially in the optical-phonon 
damping than in the dispersion law. In any case, the Frolich model has the evident 
shortcomings. 

From the experimental point of view, the best opportunity for the investigation of 
interactions between electrons and optical phonons is provided by coupled phonon- 
plasmon modes in doped semiconductors (see, e.g. ^D]). Two such modes have 
been observed in Raman experiments for many semiconductors. On the early step, 
the Raman results were compared with the theory [Tj] based on a Drude model 
(see, e.g. [12])) but more recently the Lindhard-Mermin expression for the dielectric 
function has been used [T5] . 

The Lindhard-Mermin expression ^3] represents a sophisticated generalization of 
the Lindhard function with the help of the electron relaxation time. The Lindhard 
approach is very useful while the momentum transfer k in the Raman scattering is 
compared with the Fermi momentum pp. The most essential influence of carriers 
should be expected for kvp ~ u>, where vp is the Fermi velocity and u> is the phonon 
frequency. For solids with metalic conductivity, the Fermi velocity could be estimated 
using an argument of stability under the Coulomb interaction e 2 /nTivp < 1. This 
condition gives vp ~ 0.7 x 10 8 cm/s. For the typical value of the optical phonon 
frequency uj = 500 cm -1 , the interesting values k < u/vp ~ 10 6 cm^ 1 . Therefore, the 
condition k < pp is satisfied for the carrier concentration larger than 3xl0 17 cm -3 . In 
experiments, the heavily doped semiconductors with large carrier concentration are 
used, in order to obtain a visible effect of carriers. In such a situation, the condition 
k <C pp is fulfilled, and we can apply the Boltzmann equation in calculations of 
electronic susceptibility, as well as for the evaluation of the Raman cross section. 
The method of the Boltzmann equation is valid for the anisotropic electron plasma 
in solids at arbitrary temperatures. In the present paper, we obtain the Raman 
efficiency applying the Boltzmann equation for generate carriers in heavily doped 
semiconductors when temperature is lower than the Fermi energy, T <C ef. 



2 Effective Hamiltonian and light scattering 

For the electron-phonon system in solids, we use the operator of particle numbers 
n, the phonon displacements bj, and the macroscopic electric field E, accompanying 
vibrations in polar semiconductors and acting on charges of electrons and ions. The 
effective Hamiltonian describing the inelastic light scattering in solids can be written 
in the semiclassical Wigner representation 

2 

ft = Ar / d 3 r^(r,t)U(r,t), (1) 
mc z J 

where 

M(T,t) = Tn(?,t)+g j b j (T,t)+g E E(r,t) (2) 



2 



is a linear form of variables n, bj, and E. The subscript j denotes the various 
phonon modes: longitudinal (LO) or transverse (TO). More precisely, the subscript 
j numerates the different phonon representations which can be degenerate. The 
transformation properties of the coupling constants gj are determined by this repre- 
sentation. The notation U (r, t) is introduced for a product of the vector-potentials 
of incident and scattered photons: 



A®(r,t)A w (r,t) = U(r,t) = exp[i(k ■ r - ut)]U(k, 



where the momentum and frequency transfers k = kW — k^ s \ oj = — u^ a \ The 
polarization vectors of bj(r, t), E(r, £), AW(r,t), and A^(r,t) are included in the 
coupling constants. 

The first term on the right-hand side of Eq. (J2J) descibes the light scattering by 
electron-hole pairs with a vertex 



7(P) = e^ef 



c 1 \^ ( PfnPnf PfnPnf 

<W + — —jz^ — : f „\ ; , , (i ) + 



m 



e/(p) - e n (p) + wW C/ (p) - e n (p) - uj(°) 



where the resonant term is included; and ei are the polarizatin vectors of incident 
AW(r,f) and scattered A' s '(r,t) photons. The quantum-mechanical and statistical 
average of the first term in Eq. (JSJ) 

«7^(r,t)» = / ^7(P)/ P (r,f) (3) 

can be expressed in terms of the electron distribution function f p (r, t). The constants 
gj and gs are the deformation-optic and electro-optic couplings with the phonon 
displacements and the macroscopic electric field, correspondingly. The estimation 
gives gj ~ /a 4 , gg ~ 1/ea, and 7(p) ~ m*/m, where a is the lattice parameter and 
m* is the effective mass. 

The variable U(r, t) can be considered as an external force. Then one introduces 
the generalized susceptibility x(k, oj) as the linear response to this force: 

((jV(k,oj))) =- X (k,u)U{k,u). (4) 

According to the fluctuation-dissipation theorem the function 

2 

K(k,u) = - — Imx(k,u;) 

1 — exp(— ijJjT) 

is the Fourier component of the correlation function 

K(v,t;v',t') = {{jV\v,t)M(r',t'))) (5) 
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K(k,oo)\U(k,u;)\ 2 , (6) 



which depends only on the differences r — r' and t — t' . The Raman cross section is 
given by 

r—r = — t-t I K(k. D)\U(\C. D)r. 

dwWdfiW 7rc \chmu)( l \ 

where /4 s ) is the normal to the sample surface component of the scattered wave vector 
in vacuum. 

One note should be made. Of course, any sample has the surface. The surface 
effects in the Raman scattering was considered in our work ^Hj and they are omitted 
in the derivation of Eg. (JSJ). Furthermore, the incident and scattered fields do not 
penetrate into the bulk due to the skin effect. For the optical range of the incident 
light, we have the normal skin-effect conditions. Then we integrate in Eq. (JHJ) the 
distribution \U(k, u>)\ 2 over the normal component k z . As shown in the paper the 
integration of |?7(k, u)\ 2 gives a factor 1/C2, where £2 is expressed in terms of the wave- 
vector components inside semiconductor. (2 = lm(k^' + ky). The obtained Raman 
cross section (JHJ) is dimensionless. It represents a ratio of the inelastic scattered light 
energy to the incident energy. 



3 Boltzmann equation for carriers 

The problem of the evaluation of the Raman cross section consists in the calcula- 
tion of the generalized susceptibility Q. We apply the Boltzmann equation for the 
electron distribution function: 

& J^A + v & hM. + p^M = _I [/p(r , t) _ (/pM))] . (7) 

The angular brackets denote the average over the Fermi surface 

_ 1 r, , 2dS P 



u J v(2tt) 3 ' 

where the integral is performed in the momentum space over the Fermi surface and 
u is the density of electron states, defined by the condition (1) = 1. We use the 
r— approximation which is correct for the electron scattering by defects in metals and 
heavily doped semiconductors. The collision integral in the form conserves the 
number of electrons in collisions. Therefore, the charge density satisfies the equation 
of continuity. This ensures the correct uj— dependance of the dielectric function at 
low frequencies. 

According to Eqs. (JTJ), (j2j) and instead of the unperturbed electron spectrum 
£q(p)> we introduce the local electron spectrum in the presence of the external force 
U(r, t) as follows: 

e(p,r,t) = £ (p) +7(p)tf(M) + Ci(p)&iM), 
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where the last term represents the electron-optical-phonon deformation potential 
and bj(r,t) = ((bj(r,t))). We use this form of the el-ph interaction instead of the 
Frolich poirization type C(p)divb(r, t), because it is larger by the parameter 1/ka 
for the optical phonons. 

Let us linearize Eq. ((7j), looking for its solution in the form 

f p (r,t) = / [e(p,r,t) - fi] - fr^M), (8) 

where /o[^(p ; r > t) — fi] is the Fermi-Dirac local distribution function. It is important 
that the collision term in the Boltzmann equation is canceled by the local-equilibrium 
term of Eq. (JHJ). 

We set the number-conservation condition on the chemical potential: 

/ Wf fo[£ ^ r > ') - d = / 4^ fo{£ ° - 

We obtain 

H = Ho+ (j{p))U{r,t) + (C i (p)>6 i (r,t). 
The condition implies the renormalization of vertices 

7(P)-7(P)-(7(P)), C,(P)-0(P)-(0(P)) (9) 

and this substitution should be made in that follows. 

The linearized Boltzmann equation in the Fourier components has the form 

-i(uj - k ■ v + i/r)6f p (k,uj) = ip p (k,u) + (6fp(k,u))/r, 

where 

il>p(k,u) = ev ■ E(k,w) - iu[y(p)U(k,u) + (j-(p)6j(k,w)]. 
The solution to the equation is esily obtained: 

6f p {k,u) = ity p (k,u>) + (6f p (k,u))/r]/A p , (10) 

where we denote A p = u — k • v + z/r. Now we get 

(*/ p (k,w)) =i& p (k,uj)/A p )/(l-i(T- 1 /A p )). (11) 

Notice, that in accordance with the adiabatic concept, no additional contribution 
comes from the local equilibrium distribution function fo[e(p, r, £)] in Eq. (jHJ). 
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4 Equation of motion for phonons interacting with 
carriers 



In the long- wave approximation (fc C 1/a, a being the lattice parameter), we write 
the equation of motion for the phonon displacement field 

W-^Mk,*) = ^fc«> - - ^/^(pW"). < 12 > 

where N is the number of unit cells in 1 cm 3 , M' is the reduced mass of the unit 
cell, and Z is the effective ionic charge. The nonperturbed phonon frequency ujk 
should be considered in the absence of the electric field and without any electron- 
phonon interactions. In the long-wave limit, we can expand it as lo\ = uj 2 ±s 2 k 2 with 
the value of the dispersion parameter s on the order of the typical sound velocity 
in solids. Notice that the optical phonons always have the so-called natural width 



r nat ~ ujQ^jm/M . The natural width results from decay processes into two or more 
acoustic and optical phonons. In the final expressions, we will substitute u\ — uo 2 — > 
u 2 k - iojV aat - to 2 . 

The equation (|12|) is applied to both the longitudinal and transverse phonons. It 
is seen from the Maxwell equations, that the electric field is longitudinal E || k in 
the optical region k ^> uo/c. If the excited phonons propagate in the symmetrical 
direction, the TO and LO phonons are separated. Therefore, the electric field stands 
only in Eq. (fT2j) for the LO phonon. Beside of that, the coupling Cj(p) depends on 
the phonon representation j. 

Using the solution (JTUJ), we rewrite Eq. (JTJJ): 



(OJ -J%0t,a) - |^(k, W ) = - ¥ ^ 1 , (13) 
where the phonon frequency 

uu / fflp), i«,-(p)/A p ) 2 \ 

^'-^ + MW1 V < ^7 )+ r-(VA p ) J' (14) 

the effective ionic charge 

^ v~^r ) r - (,/a p ) ) ' (15) 

and the deformation-optic coupling 

* = 9 > + uv > \ { ~~k^- ) + — ^w\) — J (16) 

are renormalized due to the electron-phonon interaction Cj(p)- 
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5 Poisson equation for the macroscopic field 



Let us consider the longitudinal electric induction D that accompanies the lattice 
vibrations. There are several contributions in the field : (1) the polarization aE(r, t) 
of the filled electron bands, (2) the lattice polarization NZb^o{ r ,t), (3) the contri- 
bution of free carrier density p = — divP e , and (4) the term P = —&HjdE = —QeU 
explicitly results from the Hamiltonian (J2J). Collecting all these terms into the 
Poisson equation, divD = 0, we find 

£oo E(k, u) + AnNZb LO (k, ^) + ^J 7^<^( k > ~ ^9 E U(k, cu) = 0, (17) 

where the high-frequency permittivity Eoo = 1 + Aira. Using the solution of the 
Boltzmann equation we rewrite the Poisson equation in the form 

e e (k, u)E(k, u) + 47rNZb LO (k, w) = 47rg E U{k, u), (18) 

where the electronic dielectric function 



e e (k,a>) = Boo + £oo^| 



kl r (u/A p (k)) 

1 - (i/AJk))/r 



with the Thomas- Fermi parameter k$ = 47re 2 i>o/e 00 . 

Due to the electron-phonon interactions Clo(p) ; the ionic charge obtains an ad- 
ditional term 

y _ y , ie "o (, Vz(lo(p) \ , i(v z /A p )(t LO (i>)/A p ) \ 

+ ~N~ [ i ~A P ~ ) ^W\) ) ' (20) 



of the opposite sign in the comparison with that in Eq. (jTojl . The electro-optic 
coupling in Eq. lJP7jl also changes, but because of the light scattering by carriers 

7(P): 

9E = 9E ~ WUo + r-(»/A p ) J • (21) 

6 Raman scattering by electron— hole pairs, phonons, 
and coupled modes 

Now we are in a position to calculate the susceptibility (@J). Using Eqs. ©, (jIUJ). 
(HH, (HD, and we get 

({M{k, u))) = - Xe (k, cu)U(k, u) + gfoik, to) + g E E(k, u), (22) 

where 

7 2 (p) v , *( 7 (p)/A p > 2 \ 



Xe(k,w) = ~UU (—7 ) + 



A p r-(i/A p ) 



gives the light scattering with excitation of electron-hole pairs. Notice, that here the 
renormalized coupling gg differs from c/e (EH) in the sign of the second term: 

,. / >*7(p) v , »(i; z /A p )(7(p)/A p ) ^ 
9e = 9e + leva — r + ttt-t-^ • 

In order to find E(k,u) and bj(\c,u), we have to solve the system of the al- 
gebraic equations (fT3J) and (fTHjh Then, using Eq. (J22|) we obtain the generalized 
susceptibility 

Y(k,a;j = Y e (k,a;j H — — . 

V ; V ; (cu]-a; 2 )£ e (k,^) + 47rArZZ/M' 

(23) 

The expression (|2*B*|) is our main result. Here the poles of the second term give 
the spectrum of collective excitations of the electron-phonon system. We discuss Eq. 
in various limiting cases. 



The electronic scattering 

We obtain the Raman electronic scattering from Eq. (f2*3|) . if we set gj = = Z = 
Z = 0. We get 

X (k,c) = Xe (k,c) + ^i T , (24) 

where §e is given in Eq. (|2*T|) with gs = 0. 

For the isotropic Fermi surface, we calculate the dielectric function Eq. (JTHJ) 
performing the integration: 

(1/A p (fc)> = In k = kv F /(oj + it- 1 ), (25) 

where one should take the branch of Inx, which is real for positive real values of x. 

For the anisotropic Fermi surface, the calculations can be done in limiting cases. 
For | ac | ^> 1, we use the expansion for the electronic dielectric function 

e e {k,u) = e oo {l + (k /k) 2 [l + i(7iv Lu/k)(6(fi)/v)]}, (26) 

where \x = v • \a/vk and 8{x) is the Dirac delta-function. In this case, the Raman 
efficiency has a "tail" due to the Landau damping: 

Im X (k,o;) = (iru uj/k)( 1 2 (p)6(i2)/v}. (27) 

We see, that the Raman cross section vanishes for the isotropic vertex 7(p) because 
of Eq. (J0|). This is a result of the Coulomb screening. It was obtained first in the 
work ^B] for the Raman scattering in semiconductors (see, |17j). 
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In the opposite case \k\ <C 1, the first term in Eq. (|24jl gives a result 

ImXeiKu) = M1 2 (P)) { ^—T , (28) 

{UT) Z + 1 

which was found first in the work ^Sj with a help of the Green's function technique. 
The second term in Eq. (J2"4~)) reveals a plasmon pole at small values of k. The 
k— expansion of the dielectric function reads 

(up- _|_ lpo\\ \ 

where the fc-independent term represents the Drude conductivity and the electron 
plasma frequency is given by the integral over the Fermi surface uj 2 e = k 2 (v 2 ). The 
complex coefficient w = k 2 ((vf) + i(v 2 ) 2 / lot) /{lo + ir^ 1 ) 2 . For the isotropic electron 
spectrum and in collisionless regime, r -1 = 0, the coefficient w = {3/h)oj 2 v F . 
The k— expansion of g F gives 

9e = -iev k(>y(p)v 2 )/(uj + zr" 1 ) 2 , 

because g F = and the zero-order term in the k— expansion vanishes due to the time 
invariance v — > —v. Then the intensity of the plasmon peak ~ k 2 in accordance with 
the known behavior of the dynamical structure factor. 



The Raman scattering by TO phonons 



Ihe second term in Eq. (|2Hj) gives the TO-phonon scattering, if we set Z = Z = gg = 

uj^q — lo z — iuL nat 

where uto, 9to are defined in Eqs. (fT4"J) and (|TT)j) with £j(p) = Cto(p); we add the 
phonon width r nat mentioned above. 

Here we should note two points. First, the TO-resonance takes place at the 
renormalized frequency £>to- Taking the real and imaginary parts of the expression 
(114)1 . we obtain the TO-phonon shift and width due to the deformation interaction 
C(p) with carriers: 

Au k = Re (o>to ~ w fc)/ 2a, A> r = rnat ~~ Im£>To/W 

Second, because of the interaction with carriers, the coupling gxo ([Hi) has an imagi- 
nary part. Therefore the line-shape of the resonance becomes asymmetric (the Fano 
resonance) : 

Tmvfk^- 1 ^TO + [(Re^To) 2 -^ 2 ]Im^ 

Im X (k,u;) _ _ [(Re - to)2 _ u2]2 - {uT)2 ■ (31) 
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The line-shape asymmetry depends on the sign of Regxo- For instance, if Re g^o > 
0, the high-frequency wing of the resonance line drops more slowly than the low- 
frequency one. In the limiting case k > 1, we expand 

9to = 9to + (^A)(7(p)Cto(p)(-^ + 2u/kv)5({i)/v) (32) 

and for k <C 1: 

9to = 9to + - . 7 P Cto P + - ( - . _ n2 • 33 

Notice, that the electron-phonon interaction Cto(p) and the light scattering 7(p) 
by carriers renormalize jointly the coupling g^Q. The frequency renormalization 

~,2 



to, 



TO 



[see Eq. (I14|) ] results only from the electron-phonon interaction Cto(p)- The 
corresponding expressions can be obtained from Eqs. (}32|) and (}33|) by substitution 
7(p) — > Cto(p)- We see, that the TO phonons become broader and harder because 
of the interaction with carriers. 

Emphasize, that the phonon renormalizations depend on the carrier density i/ 
and the average coupling Cj(p) — (0(p))- They vanish for the isotropic Fermi surface. 
The maximum value of the relative renormalization has the order of XapFm*uj/m\uj + 
ir~ l \ at kv ~ \u) + «t _1 |, where A is the dimensionless electron-phonon coupling and 
m* is effective electron mass. 



The Raman scattering by L O- phonon— plasmon coupled modes 

In this case, the carriers interact with each other and the ion vibrations via both 
the macroscopic electric field E(r,t) and the deformation potential Clo(p)- in the 
long- wave limit k — > 0, Eqs. (|15p. (|20p. and (j21j) show no renormalization of the 
ionic charge, Z = Z = Z, and the electro-optic constant, c/e = Qe- The equation 
(|2HJ) takes the form: 

goejp, u)/NM' - 4Trg 2 E (Co 2 - iuT^ - u 2 ) - ^g E ~9LoZ/M' 
X[,UJ) Xe ^> u)+ {Cj 2 -iuoT^-^)e e {^uj)+A'KNZ^M' 

(34) 

where the first term is given in Eq. (|28|L The deformation potential Clo(p) renor- 
malizes the phonon frequency Qq (fTi|. as well as the deformation-optic constant (?lo 
(|lb)). The corresponding expansion in the limiting cases are similar to Eqs. ()32|) and 
(|3*3*|) . All mentioned above about the TO line asymmetry concerns also the LO line. 
Because the dielectric function of the electron-ion system reads 

e(0, uj) = e e (0, u) + AtiNZ 2 /M\CjI - iuT nat - oo 2 ), (35) 

the second term on the right-hand side in Eq. (}34"|) has poles while e(0, uj) = 0. This 
condition defines the frequency of coupled phonon-plasmon modes in the long-wave 
limit. 
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In the absence of the electron and phonon collisions (r" 1 = r nat = 0), and 
without the el-ph interaction (C(p) — 0), one obtains with the help of Eq. the 
biquadratic equation. It gives the frequencies of the coupled phonon-plasmon modes 
at k = 0: 



I, o o , . 1 



1/2 



W ± = g He + ^LO) ± g [He + ^LO) ~ ^pe^ToJ , ( 36 ) 

where c^to = <^fc is the TO-mode frequency at k — 0, = ^to + ^pi* an d 
cjpj = 4tiNZ 2 /e^M' . These frequencies (related to the c<>ro) are shown in Fig. 
1 as functions of the electron concentration, namely, uj pe /uj^o. The upper line begins 
at u LO and tends to the electron plasma frequency u pe . The lower frequency starts 
as uj pe uj^o/uj^o and then approaches c^to- In other words, observing in the optic 
range the longitudinal phonon mode and adding electrons, we see a transition of the 
longitudinal phonon frequency from co> L o to a>ro This is a result of the Coulomb 
screening. 

We can compare Eq. (jSljl with the theory of Hon and Faust jllj . Since in that 
theory the electron-phonon interaction Clo(p) as well as the electronic scattering 
7(p) were ignored, the phonon frequency and the deformation-optic constant were 
not renormalized. Then Eq. (JMj) can be rewritten as 



X (0,u) = (4?r fff ) \ 1^(0,uj)A 2 x i /e oo - Eco/Air - 2A X i], (37) 

EooE{U, id) 

where X i = NZ 2 /M\u\ Q - iuT nat - uj 2 ), A = Cw^M'e M / 'AnN Z 2 , and C = 
gi^oZj gEMuj\ Q is the Faust-Henry coefficient. Now we see, that the expression (p?7j) 
coincides with the result of Hon and Faust [see, e.g. the paper [T3|, Eq. (3.1) ]. 

For k ytz o, Eq. (J23j) includes the dielectric function (JT5j) which differs from the 
Lindhard-Mermin expression. Second, the condition determining frequencies and 
damping of the phonon-plasmon coupled modes 

~2 • ,-pnat , .2 



(CO 



J 



iuT™ 1 - u 2 )E e (k, uj) + 4nNZZ/M' = (3£ 



contains the phonon frequency uj and the ionic charge renormalized by the electron- 
phonon interaction Clo(p)- Third, the electro-optic coupling qe ff^Tf) is modified due 
to the light scattering 7(p) from the electron-hole pairs. This effect is not canceled in 
the product QeQe of Eq. even in the absence of the electron-phonon interaction 
Clo(p)- The expansion oIqe has the form 

ieu k{v^(p)) 
9e - 9e / — . . _ n2 {ov 

for \k\ <C 1 and 

g E = 9E + eu k- 2 (uj + zr- 1 )( 7 (p)(-7r/2 - ioj/kv)5(fi)/v) (40) 
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for \k\ ^> 1. We note, that the term gp(jE has the largest imaginary part for ujt ~ 1 
and then results more in the line-shape asymmetry. 

Schematicly, the dispersion of phonon-plasmon modes is shown in Fig. 2. There 
are two main specularity of this figure. First, the behavior of the upper mode near 
the line io = kvp. Around this line (r -1 < to — kvp « kvp) the dielectric function 
(fT9j) has a singularity 



£ e (k,u) 



;.2 

'k 2 



OJ 



2kv, 



Ak 2 v 2 F 



- In — 

2 (oj - kv F ) 2 + t- 



OJ 



kvt 



(41) 



Because of this singularity, the upper mode approaches the asymptote oj = kvp, 
while the wave vector k increases. Second, in the region kvp 3> to, there is always 
one mode which mainly has phonon character. The reason of that is the decrease 
with k of the imaginary part of the dielectric function (|2fij) . 



7 Discussion 

Now we consider obtained results in a simplest way. Let us assume that the electronic 
scattering is negligibly small, 7(p) = 0. The second term in the parentheses of Eq. 
(114)1 is less than the first one in both limiting cases k < 1 and k ^> 1. We neglect 
this term at all. We also do not take into account the ionic charge renormalization 
because it vanishes at small values of k. Then we can use the expression (125)1 not only 
for the dielectric function £ e (k, oj), Eq. (|19jl . but also for the renormalized phonon 
frequency to, Eq. (JHJ). 

Results of the numerical calculations of the Raman spectra, Eq. (J2BJ), in this 
approximation are shown in Figs. 3 and 4 for two values of the electron-phonon 
coupling \ eph = z/ (£ 2 (p)) /ujqM' N ~ p F am*/m. We take values of the Faust-Henry 
coefficient C = —0.5 and the phonon natural width r nat = 10 _2 o;to- The electron 
collision rate is taken as r _1 = ICT^to, which is the usual value for heavily doped 
semiconductors JH], JH]- The wave vector k and the Thomas- Fermi parameter ko are 
given in units of cuto/^f, and the frequency oj in units of c^to i n ah figures. Both 
these figures correspond to the case of small carrier numbers u pe < c^to (see, the 
right panel of Fig. 2; for the quadratic electron spectrum oj pe = k^Vp/yZ). The left 
peak mainly has a plasmon character and the right peak is mainly the LO phonon. 
We put the ion plasmon frequency to P i = uj^o, therefore uj^o = v/Scuto- When the 
wave vector k approaches the boundary of the Landau damping region kvp > u>, 
the plasmon peak becomes broader and almost disappears at k = 0.8. The broad 
continuum in the region kvp > io results from the excitation of electron-hole pairs. 
The intensity of the plasmon peak becomes larger in the comparison to the phonon 
peak, when the electron-phonon interaction \ ev h increases. 

The A;— dispersion of the plasmon (the peak position of the Raman spectra as 
a function of k), the line width (the full width at half maximum), and the line 
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asymmetry (the frequency difference between the right and left wings of the resonance 
line at half maximum) are shown in Figs. 5-6, all in units of w rm To- The width 
and asymmetry become much larger while the plasmon peak is immersed in the 
electron-hole continuum. The maximum in this region of spectra is nothing but the 
electron-hole contribution mainly. In Fig. 6, we see for X ep h = 0.1 how close this 
maximum is located to the line to = kvF- 

The behavior of the phonon peak around uj = c^lo with increasing of k is shown 
in Fig. 7. While the wave vector increases from k — to k — 1.7, the phonon peak is 
evidently shifted to the higher frequency and becomes broader. This is effect of the 
Landau damping (see Fig. 2, the right panel). But then, k > 1.75, this peak appears 
at the lower frequency, u ~ 1.4 and becomes more sharp for k > 2.2, because the 
Landau damping decreases with k [see Eq. (J2EJ)]. 

The Raman spectra for heavily doped semiconductors and metals are shown in 
Figs. 8-9 (see also the left panel of Fig. 2). The phonon peak is located now 
around u ~ co> TO instead of uj ~ culo This is effect of the Coulomb screening: 
carriers decrease the frequency of the LO mode from cu LO to uj to . We see also that 
the el-ph interaction suppresses the phonon peak. 

The effect of the Coulomb screening and the el-ph interaction on the phonon mode 
is clearly seen in Figs. 10 - 11, where the phonon part of spectra is shown in detail. 
The lines are very asymmetric. The phonon dispersion, the line width, and the line 
asymmetry as k— functions are shown in Fig. 12. We see a singularity at k ~ ujjvp. 
It is interesting to estimate the value of the phonon dispersion. With the help of 
Fig. 11, we wind dui/dk < lO" 1 ^. From the other side, using Eqs. (jSEJ)? (f2~H|) . we 
find for the phonon dispersion uj 2 = uo\ Q + u 2 i k 2 /k 2 , which well corresponds with the 
previous estimation for our values of ko and uj p i. Note, that these estimations confirm 

the adiabatic approximation, since the value of dispersion s = uj pi /k ~ VF\Jm/M 
contains the adiabatic parameter. 

In Fig. 13, the dispersion, the line width, and the line asymmetry are shown 
for the plasmon peak in heavily doped semiconductors. Here the effect of the el- 
ph interaction on the phonon dispersion is weak, no influence on the width and 
asymmetry of the line is seen. 

8 Conclusions 

In conclusions, we emphasize first, that the result of the paper ()23|) describes the 
renormalization of the phonon frequencies, the effective ionic charge and the cou- 
pling constants due to the electron-phonon deformation interaction Cj(p)- Second, 
it involves the k— dependent semiclassical dielectric function instead of the Lindhard- 
Mermin expression. At last, the light scattering 7(p) with excitations of the electron- 
hole pairs gives not only an additional contribution Xe in Eq. (J2HJ), but also modifies 
the electro-optic qe and deformation-optic gj coupling constants which become de- 
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pendent on the frequency and momentum transfers. 
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Figure 1: Frequencies (in units of u;to) of the phonon-plasmon modes at k = in 
dependence of the free carrier concentration, namely, of the electron plasma frequency 
(in units of u TO ). We set the ion plasma frequency u pi = a>ro i n the absence of the 
free carriers. Then culo/^to — V%- 
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Figure 2: Schematic drawing of the dispersion of phono n-plasmon modes for the 
metallic {uj pe > cj to , left panel) and semiconducting {oj pe < u^o, right panel) carrier 
concentration. The dashed straight lines separate the domain kvF > uj, where the 
Landau damping exists; the dashed curves represent damped modes there. 
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Figure 3: Raman spectra from a semiconductor with low carrier concentration as 
a function of the frequency transfer uj for denoted values of the momentum transfer 
k, the Thomas- Fermi parameter k (in units of u^o/vf), and the electron-phonon 
coupling constant X ep h- We set the ion plasma frequency u pi = wto> the phonon 
natural width r nat /u; T o = 1CT 2 , and the carrier relaxation rate r _1 /r nat = 10. 
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Figure 4: Same as Fig. 3 for X eph = 0.1. 
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Figure 5: The plasmon dispersion as a function of k in units of u^o/vp (the position 
in units of o;to of the line peak of Raman spectra; upper part of the figure, solid line) . 
In the bottom, the line width (the full width at half maximum, dashed line) , and the 
line asymmetry (the difference between the right and left wings at half maximum, 
dash-dotted line) in units of cjto- The Landau damping exists to the right of the 
dotted line cu = kvp. 
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Figure 6: Same as Fig. 5 for the electron-phonon coupling \ eph = 0.1. 



21 



1 

k o = 1 


i i 






A . = 0.1 




eph 






k = 0.1 




k = 1.5 




\ / \ k = 1.7 


i 1 — 


~~\~A~— \- k-2.2 

I I 



1.3 1.4 1.5 1.6 1.7 1.8 1.9 

frequency (rel. un.) 



Figure 7: The LO phonon Raman spectra for large momentum transfers k. 
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Figure 8: Raman spectra from a heavily doped semiconductor; notations as in Fig. 
I. 
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Figure 9: Same as in Fig. 8 for X eph = 0.2. 
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Figure 10: The LO phonon part of the Raman spectra from heavily doped semicon- 
ductor for various momentum transfers. 
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Figure 11: Same as in Fig 10, but the el-ph interaction is taken into account. 
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Figure 12: Dispersion of the phonon peak (upper), the line width and the line asym- 
metry (bottom). The boundary of the Landau damping region is shown as a dotted 
line. 
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Figure 13: Dispersion of the plasmon peak. 
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